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CO " ABSTRACT 

<N 

Context. Conceptually reconstructing the physical conditions in relativistic jets, given the observed electromagnetic spectrum, is a 
complex inverse problem. 

Aims. We aim to improve our understanding of the mechanisms operating in relativistic jets by modeling their broadband electro- 
magnetic spectrum. 

Methods. We develop an inhomogeneous jet model where the injection of relativistic primary and secondary particles takes place in 
i-C " a spatially extended region. We calculate the contribution of all particles species to the jet emissivity by several radiative processes, 

and assess the effect of gamma-ray absorption in internal and external photon fields. A number of specific models with different 
parameters are computed to explore the possibilities of this scenario. 

Results. We obtain a variety of spectral shapes depending on the model parameters, some of them predicting significant gamma-ray 
emission. The observed broadband spectrum of the low-mass microquasar XTE Jl 118+480 can be satisfactorily reproduced by the 
model. 

Conclusions. Our results indicate that outbursts similar to those displayed in the past by XTE Jl 118+480 might be detected with 
present-day gamma-ray instruments. 

> 
O 

WV 1. Introduction the Fe Ka line at ~ 7 keV in some XRBs (e.g. Miller et al. 

2002, Tomsick et al. 2009, Duro et al. 2011). The transition 
^ . Jets are collimated outflows. Astrophysical objects that launch from an outer thin disk to an inner quasi-spherical corona is 
t-H ; jets have two common components: a rotating compact object se lf-consistently predicted in some accretion models, such as 
! and an accretion disk. Examples of these systems are gamma- t h e advection-dominated accretion flows (ADAFs) studied by 

■ ray bursts, active galactic nuclei, microquasars, and young stellar Narayan & Yi ( 1 994, 1 995) and the adiabatic inflow-outflow so- 
^ ! objects. lutions (ADIOS) of Blandford & Begelman (1999). 

• i-h ■ Microquasars are X-ray binaries (XRBs) that produce rela- 
" tivistic jets. The jets of microquasars have bulk Lorentz factors 
H . of ~ 1 — 10, and a typical power of 10 35-37 erg s _1 (e.g. Gallo 

■ et al. 2005, Sell et al. 2010, see also Heinz & Grimm 2005 and 
references therein). Steady jets in microquasars appear in the so- 
called low-hard spectral state. In this state, the X-ray spectrum is The ejection of jets during the low-hard state is revealed 
the sum of a low-luminosity black body component and a hard by the flat radio spectrum, which is interpreted as synchrotron 
power-law that cuts off at ~ 100 - 150 keV. These characteris- radiation of non-thermal electrons in an expanding outflow 
tics are well reproduced by "disk+corona" models (see for ex- (Blandford & Konigl 1979). Broadband observations of some 
ample Poutanen et al. 1997, Dove et al. 1997, and Malzac et al. XRBs indicate that there is a correlation of the form 5 ra di x S 1 ^ 7 
2001). In these models, the thermal emission originates in an ac- between the radio and the X-ray flux (Gallo et al. 2003). This 
cretion disk (truncated at tens or hundreds of gravitational radii suggests that the jets might contribute significantly to, or even 
from the compact object), and the hard power-law is generated dominate, the X-ray emission in the low-hard state. Two con- 
in a "corona" of hot plasma that surrounds the black hole. The firmed microquasar - Cygnus X-l (McConnell et al. 2002, Albert 
main radiative mechanism of the corona is the Comptonization et al. 2007, Sabatini et al. 2010, Laurent et al. 201 1) and Cygnus 
of disk photons by thermal electrons or electron-positron pairs. X-3 (Abdo et al. 2009, Tavani et al. 2009) - have been detected 
The presence of the corona can also explain the detection of in gamma rays. The jets - loaded with high-energy particles and 
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moving at relativistic bulk velocities - appear as the probable site 
of the gamma-ray emissionQ 

Jet radiative models have been typically divided into "lep- 
tonic" (e.g. Markoff et al. 2001, Kaufman Bernado et al. 2002, 
Bosch-Ramon et al. 2006) and "hadronic" (e.g. Romero et al. 
2003, 2005). In hadronic (leptonic) models, the bulk of the 
gamma-ray emission is produced in interactions initiated by rela- 
tivistic protons (electrons). Although microquasar jets can prop- 
agate over distances of hundreds of AU (Mirabel & Rodriguez 
1999), most of these models are based on the "one-zone" ap- 
proximation: the emission region is assumed to be homogeneous 
and compact. 

In this work, we present an inhomogeneous lepto-hadronic 
model for the radiation of jets. We intend to expand available 
models by adopting more general assumptions about the con- 
tents of relativistic particles and the characteristics of the acce- 
leration/emission region. In our model, relativistic protons and 
electrons are injected in a spatially extended and inhomogeneous 
region of the jet. We investigate the consequences of the prop- 
agation and cooling of these particles (and also of secondary 
muons, pions, and pairs) on the jet radiative output. The charac- 
teristics of the spectral energy distributions (SEDs) can provide 
insight into the physical conditions in the jets. As an applica- 
tion, we fit the spectrum of the low-mass XRB XTE Jl 1 18+480, 
a very well-studied black hole candidate in the Galactic halo. 

The paper is organized as follows. In Sect. [2] we introduce 
the main features of the model. We present the general assump- 
tions made about the geometry and energetics of the jet, and 
briefly describe the interaction processes between the relativistic 
particles and the matter, radiation, and magnetic field. We out- 
line the method of calculation of the particle distributions and 
the jet radiative spectrum, including absorption corrections. In 
Sect. [3] we exhibit and analyze the results of a series of models 
calculated for different sets of parameters. Section|4]presents the 
application to XTE Jl 1 18+480. We show the best-fit SEDs and 
discuss the possible detection of the source with present and fu- 
ture gamma-ray telescopes. We close with a short review of the 
model results and the perspectives for future studies. 




low-mass 
donor star 



Fig. 1. Sketch of a low-mass microquasar. 
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Fig. 2. Detail of the jet and the acceleration region (not to scale). 
Some relevant parameters are indicated. 



2. Inhomogeneous jet model 

2.1. The jet 

We generalize the one-zone model presented in Romero & Vila 
(2008), Reynoso & Romero (2009), and Vila & Romero (2010) 
to treat the injection and cooling of relativistic particles in an 
extended inhomogeneous region of the jet. 

The main components of the binary system are sketched in 
Fig. [T] a sketch of the jet is shown in Fig. [2] A compact object 
(hereafter a black hole) of mass Mbh accretes matter from a low- 
mass companion star that has filled its Roche lobeQ We write the 
accretion power L accr in terms of the Eddington luminosity of the 
black hole as 



iaccr = Mc - q aca - L E dd 



1.3 X 10 38 4 a 



Mi 



M, 



BH 1 -i 

erg s , 



(1) 



1 If a fraction of the particles in the corona are relativistic, the corona 
might also radiate gamma rays; see Romero et al. (2010) and Vieyro & 
Romero (in preparation). 

2 In low-mass microquasars, the donor star is old and dim, and the 
stellar wind and radiation field are negligible targets for the relativistic 
particles in the jet. 



where M is the mass accretion rate, M is the solar mass, and 
Qaccr is an adimensional parameter. An accretion disk extends 
from an inner radius R[ n to an outer radius R oa . Nearer the 
black hole, the plasma "inflates" and forms a hot, optically thin, 
corona. 

A fraction of the accreted matter is ejected into two symmet- 
rical jets, each carrying a power 



_ 1 

^jet — ^ ^/jet ^acci'7 



(2) 



where qj et < 1 . The jets are launched perpendicularly to the plane 
of the accretion disk at a distance zo from the black hole, and 
propagate up to z = Zend- The jet axis makes an angle #j et with 
the line of sight of an observer on Earth. The outflow expands as 
a cone as it advances, hence its radius grows as 



(3) 



r(z) = r - . 



None of the parameters of the model depend on the radial coor- 
dinate (radial symmetry assumption). 

If the plasma is ejected by some magnetocentrifugal mecha- 
nism, the magnetic energy density at the base of the jet must be 
enough to set the plasma in motion. We fix the value Bq of the 
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magnetic field at zo by equipartition between the magnetic and 
the kinetic energy densities 



t^mag(Zo) = £4in(Zo)- 



(4) 



The magnetic energy density decreases as the jet expands. We 
parameterize the dependence on z of the magnetic field strength 

as 



(5) 



where m > 1 (e.g. Rrolik 1999). 

In the usually accepted model of jet acceleration (see for ex- 
ample Spruit 2010), a fraction of the magnetic energy is con- 
verted into kinetic energy of the plasma. The bulk Lorentz factor 
of the jet, rj et , then increases up to a terminal value. The behav- 
ior of Fjet with the distance to the black hole can be studied both 
analytically and numerically using the equations of the MHD 
(e.g. Lyubarsky 2009, Tchekhovskoy et al. 2008, 2010); a sim- 
pler approach is presented in Reynoso et al. (2011). Here we 
simply adopt a constant value of Tjet for all z, leaving dynamical 
considerations for forthcoming works. 

In the region z acc < z < z m ax, some fraction of the jet power 
is transformed into kinetic energy of relativistic electrons and 
protons. We assume that the acceleration mechanism is particle 
diffusion through shock fronts, which is also known as diffusive 
shock acceleration. This assumption constrains the value of z acc , 
since for shock waves to develop the magnetic energy density 
of the plasma must be in sub-equipartition with the kinetic en- 
ergy density (Komissarov et al. 2007). We then fix z acc from the 
condition 



^mag(^acc) — P^kin(^acc)» 



(6) 



where p < 1 . The power transferred to the relativistic particles is 

L re l = qrelLjet, (7) 

where qj et < 1. The value of L re i is the sum of the power injected 
in both protons and electrons 



In addition, we assume that L p and L e are simply related as 



(8) 



(9) 



If a > 1 the jet is proton-dominated, whereas for a < 1 most of 
the power is injected in relativistic electrons. 

2.2. Particle cooling and maximum energy 

The relativistic particles in the jet lose energy by radiation and as 
they exert adiabatic work on the walls of the expanding plasma. 
Several processes contribute to the radiative cooling, since par- 
ticles can interact with the magnetic, radiation, and matter fields 
of the jet. 

The most important interaction channel with the jet magnetic 
field is synchrotron radiation. For a charged particle of energy E 
and mass m, the synchrotron energy-loss rate is (e.g. Blumenthal 
& Gould 1970) 



dE_ 

dt 



synchr 3 \ ttl I 



4 (m e \ 3 CCT T JJ. 



mag 



n 2 ' 



(10) 



where c is the speed of light, m e is the electron mass, cr T is the 
Thomson cross-section, and U mdg is the magnetic energy den- 
sity. The ratio (m e /m) 3 in Eq. (TlOb makes synchrotron cooling 
particularly efficient for the lightest particles. 

Relativistic electrons also cool through inverse Compton 
(IC) scattering off the jet radiation field. We calculate the energy- 
loss rate in the Klein-Nishina regime following Blumenthal & 
Gould (1970). We only consider the synchrotron field of primary 
electrons as a target field. To estimate the density of synchrotron 
photons, we apply the "local approximation" of Ghisellini et al. 
(1985) 



^ synchr ( £"? Z) ~ £synchr(^ Z) 



r 3 et(z) 



(11) 



where s syncm -(e, z) is the synchrotron power per unit volume per 
unit energy at energy e. We do not consider the IC cooling of 
protons, since the energy-loss rate for this process is much lower 
than that of other channels of interactions between energetic 
hadrons with radiation. 

The interaction of relativistic protons with radiation can cre- 
ate electron-positron pairs ("photopair" production) 



p + y — > p + e + e 



(12) 



The photon threshold energy is ~ 1 MeV in the rest frame of 
the proton. At higher energies, proton-photon collisions can also 
yield pions ("photomeson" production). The two main channels 
are 



p + y — > p + got + b (n + + n~) 
and 

p+y— >h+7t + + cm + b(n + + n~) , 



(13) 



(14) 



where the integers a and b are the pion multiplicities. Both chan- 
nels have approximately the same cross-section, and a photon 
threshold energy of ~ 145 MeV in the proton rest-frame. 

The subsequent decay of charged pions injects leptons and 
neutrinos 



p. +v fl , fx 
pT +v ll , p 



+ v e + v„ 



+ V e + V„ 



whereas neutral pions decay into two photons, 



y + y. 



(15) 
(16) 

(17) 



For photomeson and photopair production, the proton 
energy-loss rate can be calculated as in Begelman et al. (1990). 
Once again, we only consider as a target radiation field the syn- 
chrotron photons of primary electrons calculated with Eq. ( fTTT ). 

High-energy protons can also interact with the thermal par- 
ticles in the jet plasma. If the energy of the relativistic proton is 
higher than the threshold value for n° production (« 1 .22 GeV), 
the collision with a thermal proton can create pions 



p + p^>p + p + an + b (n + + n ) , 

p + p — > p + n + 7r + + an + b (n + + n~) . 



(18) 
(19) 



As in proton-photon interactions, proton-proton collisions inject 
photons by means of the decay of neutral pions. Charged pions 
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inject electron-positron pairs and neutrinos through the decay 
chains of Eqs. (TT5b and ( fT6b . 

The energy-loss rate for a proton of energy E is given by 
(e.g. Begelman et al. 1990) 



dE 
dt 



Yip C E KppCTpp, 



(20) 



/'/' 



where n p is the number density of thermal protons, o- pp is the 
proton-proton cross-section, and k pp a 0.5 is the total inelastic- 
ity of the interaction. A convenient parametrization for cr pp is 
presented in Kelner et al. (2006). We calculate the value of n p in 
the observer rest-frame^ as in Bosch-Ramon et al. (2006) 



n p = (l - <2jet) : 



-jet 



' jet 



2 ' 



(21) 



where m p is the proton mass. 

Leptons can cool efficiently by interacting with the ther- 
mal protons in the jet through relativistic Bremsstrahlung. The 
energy-loss rate for a particle of energy E, charge e, and mass m 
is 



dE 
dt 



= -4 r] a FS c n p (z) 





r 




MB) 




E, 


\ mc L j 


~ 3 



(22) 



the physical conditions in the jet are such that the cooling time 
is shorter than the mean lifetime. 

Our treatment of the cooling of secondary particles includes 
adiabatic and synchrotron losses. For muons and electron- 
positron pairs, we also consider IC scattering, and pion-photon 
and pion-proton collisions for charged pions. The cross-section 
of pion-proton interactions is w 2/3 of that for proton-proton 
collisions (Gaisser 1990). We apply the same approximation for 
pion-photon collisions. 

The injection of secondary particles is not a consequence of 
diffusive shock acceleration, but a by-product of the interaction 
of primary protons and electrons. The maximum energy of pi- 
ons, muons, and secondary electron-positron pairs is then not 
fixed by Eq. (l24l i. It is instead determined by the characteristics 
of the particular process through which they were injected, and 
by the maximum energy of the primary particles. 

2.3. The steady-state distribution of relativistic particles 

We parameterize the particle injection function of primary elec- 
trons and protons in the jet reference frame (in units of erg 1 
s cm 4 ) as a power-law in energy multiplied by an exponential 
cutoff 



Q(E,z) = Q E- a exp[-E/E maiL (z)] f(z). 



(27) 



where r e is the classical electron radius and o-ps is the fine struc- 
ture constant. 

Finally, to calculate the adiabatic energy loss rate we use the 
expression given in Bosch-Ramon et al. (2006) 



dE 
dt 



(E,z) 



ad 



2 Vjet 

"3T J 



(23) 



The maximum energy of primary electrons and protons at z can 
be estimated by equating the total energy-loss rate to the accele- 
ration rate 



dE 
dt 



t dt 



max ■> z) - 



(24) 



The total energy-loss rate is simply the sum of the relevant loss 
rates for the particle species in consideration 



dE 
dt 



(E,Z) = Y J 



cTE 
dt 



(E,z). 



(25) 



The acceleration rate is (e.g. Aharonian 2004) 



dE_ 

dt 



= rjecB(z), 



(26) 



where e is the electron charge, and r\ < 1 is a phenomenologi- 
cal parameter that characterizes the efficiency of the acceleration 
mechanism. 

Pions and muons can also interact before decaying. This ef- 
fect becomes important when particles are very energetic and 



3 Henceforth, we refer to as "comoving" or "jet frame" the reference 
frame fixed to the jet plasma, where the outflow bulk velocity is zero. 
Its z-axis is parallel to the jet symmetry axis. The "observer frame" 
correspond to that where the jet bulk velocity is Vj ct . When it is necessary 
to distinguish between the same quantity in both frames, the value in the 
jet frame is denoted by primed symbols. 



For diffusive shock acceleration, the value of the spectral index 
is typically in the range 1.5 < a < 2.4 (Rieger et al. 2007). The 
cutoff energy E max depends on z, and is calculated by balancing 
the energy-loss rate and the acceleration rate, as explained in 
Sect. 12.21 The function f(z) is a step-like function that tapers the 
length of the acceleration region 



m = i 



i 



1 + exp[-(z - z max )] 



1 Z ^ Zmax? 
Z > Z max - 



(28) 



Finally, the value of the normalization constant Qo is obtained 
from the total power injected in each particle species 



U = f dr f °° 

JV JE mi „ 



dE E Q(E, z), 



(29) 



where the index i runs over protons and electrons. 

As shown in Sect. 12.21 the interaction of relativistic protons 
with matter and radiation injects charged pions, muons, and sec- 
ondary electron-positron pairs. Electron-positron pairs are also 
created through photon-photon annihilation. The expression of 
Q(E, z) for secondary particles depends on the specific injection 
process. 

The sources of charged pions are proton-photon and proton- 
proton interactions. Suitable expressions for the corresponding 
pion injection functions were presented in Atoyan & Dermer 
(2003), Kelner et al. (2006), and Kelner & Aharonian (2008). 
We use the formulas of Lipari et al. (2007) for the muon injec- 
tion function. Finally, to calculate the electron-positron injection 
rate from muon decay we follow Schlickeiser (2002). 

To estimate the pair output from photopair production, we 
apply the formulas of Chodorowski et al. (1992) and Miicke et 
al. (2000). The last source of pair injection that we consider is 
photon-photon annihilation 



r + 7 ■ 



e + 1 



(30) 
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The pair injection function for this process is derived in Bottcher 
& Schlickeiser (1997). Two-photon annihilation is also a photon 
sink, and can eventually modify the production spectrum. Its ef- 
fect on the escape of photons from the jet is discussed later. 

Once the injection function is known, we calculate the 
isotropic, steady-state particle distributions N(E, z) (in units of 
erg -1 cm -3 ) in the jet reference frame by solving a transport 
equation of the form (e.g. Khangulyan et al. 2008) 



dN d idE 



\ N 
N\ + 

tot / Tdec(E) 



Q(E,z). 



(31) 



Tyy(Ey,Z) = ~ J J J (Tyy(Ey, €, U) «(&, f) X 

(1 — u)de du dl, (36) 

where I is the length of the path traversed by the photon from the 
emission site to the observer, n(e, r) is the energy distribution of 
targets photons, u is the cosine of the collision angle, and CJ~y-y IS 
the annihilation cross-section (e.g. Gould & Schreder 1967). The 
threshold energy of the target photon depends on the collision 
angle, and is given by 



This equation is appropriate to describe particle distributions 
over spatially extended regions. It incorporates the effects of the 
variation with z in the parameters that govern the interaction of 
the relativistic particles, namely the magnetic field, the radiation 
fields, and the density of thermal particles. From left to right, 
the three terms on the left-hand side of Eq.PTTl account for the 
changes in N{E,z) caused by particle convection along the jet, 
energy losses, and decay. The convection velocity is on the order 



of the jet bulk velocity, 



"jet 



The decay term is nonzero 



only for pions and muons, where Td ec (£') is the corresponding 
mean lifetime in the jet frame. 

2.4. Non-thermal spectrum of the jet 

We calculate the radiative spectrum of all particle species pro- 
duced by each of the interaction processes described in Sect. 
12.21 Detailed formulas can be found in Romero & Vila (2008), 
Reynoso & Romero (2009), and references therein. 

The spectra are first obtained in the co-moving frame, and 
then transformed to the observer frame. If q'y is the jet-frame 
photon emissivity (in erg s cirT 3 ) associated with process i, the 
specific luminosity at energy E' is 



l;«(£;) = 2n 



£ 



rf et (z)qf(E'z)dz. 



jet 



(32) 



e±(Ey, u) = 



Ey (1 - U) 



(37) 



Eq. (l36*i l is simplified when the target photon distribution is 
isotropic. A convenient expression for t J7 in this particular case 
is obtained in Gould & Schreder (1967). 

3. General results 

We compute several models with different values of some of 
the parameters. Here we present the results for four represen- 
tative models. In all cases, we fix M BH = 10M o , zo = 5Otf grav 



r = 0.1 zo, Zend = 10 12 cm, 6L, = 30°, q 



0.1, 



<7jet 



0.1, 



q K \ — 0.1, T] — 0.1, and E^ = 10 m^^c . The parameters spe- 
cific to each model are listed in TableQ] 



3.1. Cooling rates 

In Figs. [3] and |4] we plot the cooling rates for primary electrons 
and protons in model A. They are calculated at z = z aC c (base of 
the acceleration region), z = z m ax (top of the acceleration region), 
and z = Zend ("end" of the jet). The cooling (acceleration) rate for 
process i is defined as 



_ t _ I dE 
f/ ~~E ~dt 



(38) 



The luminosity in the observer frame is then 
L®(E 7 ) = D 2 LfXE' r ), 



where 



Ey = DE' y 



is the photon energy in the observer frame, and 

-l 



D = [r jet (l -jSjetCOSejet)] 



(33) 



(34) 



(35) 



is the Doppler factor. In Eq. (l35l l. # je t is the jet viewing angle and 



y8 



'jet 



'jet 



jc. In the particular case of proton-proton interactions, 



the luminosity is calculated directly in the observer frame. Since 
Eq. ( l3TT l yields the proton distribution in the jet frame, it must 
be first transformed to the observer frame. For this, we use the 
general expressions given in Torres & Reimer (201 1). 

The emission spectrum must be corrected for absorption 
caused by photon-photon annihilation into electron-positron 
pairs (see Eq. d30b). The optical depth for a photon of energy 
Ey created at height z on the jet axis is 



Synchrotron radiation dominates the cooling of electrons near 
the base of the acceleration region. Further away from the 
black hole, this process gradually becomes less relevant at lower 
energies, since t^ nchi oc z~ 3 . The maximum energy of elec- 
trons, however, is always determined by the balance of the syn- 
chrotron cooling rate and the acceleration rate. Synchrotron-self- 
ComptorQ (SSC) energy losses are in general much smaller, and 
the cooling caused by relativistic Bremsstrahlung is negligible. 

Adiabatic energy losses are the most important for protons 
all along the jet. We note that the cooling produced by proton- 
photon (py) interactions is completely negligible at large z. 



3.2. Particle injection and distributions 

Figure [5] shows the dependence of the injection function on en- 
ergy and z, for primary electrons and protons in model A. The 
injection is confined to the region z < z max = 10 1() cm. As ex- 
pected from the cooling rates in Figs. [3] and [4] protons reach 
energies much higher than electrons. The maximum energy of 
electrons is determined by the synchrotron losses, and hence it 



4 iJgtav is the gravitational radius of the black hole, R srm = GMbh/c 2 . 

5 Synchrotron-self-Compton (SSC) scattering is the IC scattering of 
particles's own synchrotron radiation field. 
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Table 1. Values of the model parameters for four representative models. 



Parameter (symbol) 


Model A 


Model B 


Model C 


Model D 


Magnetic field decay index (in) 


l.j 


Z.U 


L.J 


1 c 
L.J 


Ratio U B /U k atz acc (p) 


0.9 


0.9 


0.5 


0.1 


Base of the acceleration region (z acc ) 


1.8 x 10 s cm 


1.2 x 10 s cm 


3.2 x 10 8 cm 


1.6 x 10" cm 


End of the acceleration region (z max ) 


10 10 cm 


10'° cm 


10" cm 


10'° cm 


Ratio L p /L e (a) 


1 


100 


1 


1 


Power relativistic protons (L p ) 


3.2 x 10 35 erg s _I 


6.4 x 10 35 erg s _I 


3.2 x 10 35 erg s _I 


3.2 x 10 35 erg s -1 


Power relativistic electrons (L„) 


3.2 x 10 35 erg s -1 


6.4 x 10 33 erg s _1 


3.2 x 10 35 erg s"" 1 


3.2 x 10 35 erg s" 1 


Particle injection spectral index (a) 


1.5 


1.5 


1.5 


2.2 




7 8 9 10 11 12 13 14 7 8 9 10 11 12 13 14 7 8 9 10 11 12 13 14 

Log 10 (E e /eV) Log m (E e /eV) Lo S 10 ( E e /eV ) 



Fig. 3. Cooling and acceleration rates for relativistic electrons for model A at the base (left) and the top (center) of the acceleration 
region, and at the "end" of the jet (right). 
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11) p 11) p It) p 



Fig. 4. The same as in Fig.[3]but for relativistic protons. The abbreviations and "py" stand for proton-proton and proton-photon, 
respectively. 



increases with z as the magnetic field decreases. For protons, 
adiabatic losses are the main cooling channel, and the maximum 
proton energy decreases with z. 



The steady-state particle distributions calculated from Eq. 
( l3TT i are plotted in Fig. [6] We note that the most energetic elec- 
trons quickly disappear after the injection is switched off - they 
cool and accumulate at lower energies. Since the cooling times 
are much longer for protons than for electrons, the behavior of 
the proton distribution is quite different. The number of the most 
energetic protons also decreases forz > z max , but there are plenty 
of high-energy protons outside the acceleration region. In one- 
zone models, where there is no convection term in the transport 
equation, or it is replaced by an escape term, the particle distri- 
butions would be zero in the region where Q(E , z) ~ 0. 



3.3. Spectral energy distributions 

Figure [7] shows the SEDs obtained for the four models of Table 
[TJ which have a broad range of spectral shapes. 

In model A, the power injected in relativistic protons is 
large (a = 1). The jet emission up to ~ 1 TeV is synchrotron 
and IC radiation of primary electrons, reaching luminosities of 
~ 10 36 erg s~' at ~ 10 MeV The very high-energy tail of the 
spectrum is due to proton-proton and proton-photon interactions. 

In model B, most of the energy is transferred to relativis- 
tic protons (a = 100). The synchrotron radiation field of primary 
electrons is considerably weaker, and so are all other interactions 
for which this photon field is a target (proton-photon collisions 
and IC scattering). Furthermore, in this model m = 2 and the 
magnetic field strength decreases rapidly with z. This also con- 
tributes to quenching the synchrotron emissivity. We note that 
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Fig. 6. Energy distribution of relativistic electrons (left) and protons (right) for model A. Figures available in color in the electronic 
version of the manuscript. 



the radiative output of proton-proton collisions is only slightly 
affected compared to that of model A. 

The only difference between model C and model A is the 
extent of the acceleration region. In model C, the base of the re- 
gion is shifted to slightly higher z, and extends up to 10 11 cm. 
This "spread" in the spatial distribution of the relativistic parti- 
cles affects the proton-proton gamma-ray spectrum. The more 
extended the acceleration region, the less radiatively efficient 
this process is. 

Finally, in model D the injection spectral index of the rela- 
tivistic particles changes from a — 1.5 to a = 2.2. This turns 
the electron synchrotron spectrum from hard to relatively soft. 
The same happens with the gamma-ray emission produced by 
proton-proton collisions. 

Both Bremsstrahlung and synchrotron radiation of sec- 
ondary particles do not contribute significantly to the SED of 
any model. 

The total luminosity corrected by absorption is also plotted 
for every model of Fig. [7] The effect of absorption is only notice- 



able at very high energies in models A and D, where the accele- 
ration region is less extended. 

These four models do not exhaust the possible parameter 
space that can be explored. To calculate the SEDs of models 
A to D, we choose a set of fixed parameters and study the ef- 
fects of varying the rest of them. Modifying other parameters, 
such as <7 accr or 77, would introduce further interesting changes in 
the spectrum, including higher gamma-ray luminosities. This is 
clearly exemplified by the SEDs calculated to fit the spectrum of 
XTE Jl 1 18+480 presented in Sect. 1331 



4. The low-mass microquasar XTE J1 118+480: 
characteristic parameters and observations 

XTE Jl 1 18+480 is an X-ray binary in the Galactic halo. It hosts 
a low-mass donor star (M, « 0.37M o ) and a black hole of mass 
M BH ~ 8.53M (Gelino et al. 2006). 

The estimated distance to XTE Jl 118+480 is d * 1.72 kpc 
(Gelino et al. 2006). The source is located at high Galactic lati- 
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Fig. 7. Spectral energy distributions for the four general models of Table Q] The various curves labeled "synchr." correspond to the 
synchrotron radiation of primary electrons (the most luminous component), protons, and secondary particles (pions, muons, and 
electron-positron pairs). Figures are available in color in the electronic version of the manuscript. 



tude (b = +62°) in a region of low interstellar absorption along 
the line of sight. This peculiarity of its position allows to perform 
very clean observations. 

XTE Jl 118+480 is a transient XRB, spending long periods 
in quiescence. Since it was first detected in 2000 (Remillard et al. 
2000), two outbursts were observed: at the time of its discovery 
and in 2005. The source was extensively observed at different 
wavelengths during the two episodes, and on both occasions the 
spectrum showed the characteristics of the low-hard state. No 
outflows were directly imaged, but the presence of jets can be 
inferred from the radio-to-infrared/optical emission (Fender et 
al. 2001). 

The outburst of 2000 lasted for about seven months. 
Observational data with very complete spectral coverage were 
presented and analyzed in Hynes et al. (2000), McClintock et al. 
(2001), Esin et al. (2001), and Chaty et al. (2003). As shown in 
Chaty et al. (2003), the SED did not change significantly over 
a period of about three months. The second outburst started in 
January 2005 (Zurita et al. 2005; Pooley 2005; Remillard et al. 
2005) and lasted for 1-2 months (Zurita et al. 2006). Radio-to- 
X-ray data from this epoch were presented in Hynes et al. (2006) 
and Zurita et al. (2006). 



The radio-to-X-ray spectrum of XTE Jl 1 18+480 in outburst 
has been interpreted as the sum of the black body emission of a 
thin accretion disk and synchrotron radiation from non-thermal 
electrons in a jet (see, for example, the works of Markoff et al. 
2001, Yuan et al. 2005, Maitra et al. 2009, and Brocksopp et al. 
2010). Other authors replace the contribution of the jet by that of 
a hot, optically thin ADAF (see Esin et al. 2001 and Yuan et al. 
2005). The ADAF models, however, underpredict the observed 
radio and UV emission. 

We apply our inhomogeneous jet model to fit the broad- 
band data from the two known outbursts of XTE Jl 118+480. 
Differentiating from previous works, we explore the conse- 
quences of the injection of non-thermal protons and secondary 
particles in an extended region of the jet. We also consider the 
effects of internal and external absorption on the jet emission. 

Table|2]shows a brief log of the observations used for the fits. 
The data from the 2000 outburst are taken from McClintock et al. 
(2001), whereas for the 2005 outburst we use the data published 
in Maitra et al. (2009)0 Please refer to these works for details 
of the instrumental techniques and the reduction process. The 



6 We extracted the data with the help of the ADS's Dexter Data 
Extraction Applet and a script prepared by the authors. 



8 



Vila et al.: An inhomogeneous lepto-hadronic jet model 



UV/X-ray spectrum of the 2000 outburst displays a "dip" in the 
energy range 0.15-2.5 keV. According to Esin et al. (2001), this 
feature might be caused by absorption in a region of partially 
ionized gas interposed in the line of sight. Following Maitra et 
al. (2009), we exclude this energy band from the fit. 



Table 2. Log of the observational data used in the fits. 



Observation date 


Instrument 


Range 


2000 April 18 


Ryle Telescope 


15.2 GHz 




UKIRT" 


1-5 fim 




HST i 


1155-10250 A 




EUVE r 


0.1-0.17 keV 




Chandra 


0.24-7 keV 




RXTE - PCA rf 


2.5-25 keV 




RXTE - HEXTE e 


15-200 keV 


2005 January 23 


Ryle Telescope 


15.2 GHz 




UKIRT 


J, H, K-band 




Liverpool Telescope 


V, B-band 




RXTE 


3-70 keV 



Notes. <o) United Kingdom Infrared Telescope {b) Hubble Space 
Telescope (c) Extreme Ultraviolet Explorer w Rossi X-Ray Timing 
Explorer - Proportional Counter Array w Rossi X-Ray Timing Explorer 
- High Energy X-ray Timing Experiment 



5. Fits to the SED of XTE J111 8+480 in low-hard 
state 

5.1. Accretion disk model 

The shape of the observed SED of XTE Jl 1 18+480 in the IR- 
optical/UV band strongly suggests that it is produced by thermal 
radiation from an accretion disk; we use a simple model to fit 
this component. We parameterize the radial dependence of the 
temperature as (e.g. Frank et al. 2002) 



T(R) 



0.488 \R 



R 



-3/4 




(39) 



where R[ R is the disk inner radius, and r max is the maximum 
disk temperature reached at R max = (49/36)/?; n . This tempera- 
ture profile is consistent with that of a standard optically thick, 
geometrically thin accretion disk (Shakura & Sunyaev 1973). 

Every annulus of the disk radiates as a black body at the local 
temperature T(R). The observed flux at energy E y is then 



F d (E 7 ) = 2n 



cos 6>d 



B(E 7 ,R)RdR, 



(40) 



where B(E y , R) is the Planck function 



B(E 7 ,R) 



" 2 h 3 exp [E y /kT(R)] 



(41) 



and #d is the inclination angle of the disk with respect to the line 
of sight. We assume that the disk is perpendicular to the jet, so 

The power emitted per unit area of the disk is 
D(R) = <r SB T{R) 4 , where <xsb is the Stefan-Boltzmann constant. 
The total luminosity of the disk can be calculated by integrating 
D(R) over the two faces of the disk 



Li = 2 x 2ncrs\ 



T(R) R dR 



An j 7 

"Mo 



.488 



(42) 



The approximation is valid for R mX » R m - 

In steady state, half of the gravitational energy lost by the 
infalling matter is radiated in the disk. Then 



L _ 1 GM BH M _ 1 R 
2 R m 2 R n 



grav • 9 

Mc 1 . 



(43) 



Eqs. (l42b and (l43l provide an estimation of the accretion power 



2 8;r / T max \ 4 Ri 



(44) 



The photons from the disk irradiate the jet. Since the disk is 
considerably luminous, the interaction of the disk photons with 
the relativistic particles in the jet cannot be neglected. We con- 
sider the IC scattering of the disk photons by the primary elec- 
trons and muons in the jet. In the observer frame, the interaction 
is anisotropic, and the dependence of the Klein-Nishina cross- 
section on the collision angle must be considered. We calculate 
the anisotropic IC spectrum as in Dermer & Schlickeiser (1993). 
To simplify the formulas, the disk radiation field at fixed radius is 
approximated as a monoenergetic field at energy (e) = 2.7kT(R), 
where k is the Boltzmann constant. Under the same approxima- 
tions, an analytical expression for the energy-loss rate for this 
process is obtained in Bottcher et al. (1997)0 

High-energy photons emitted in the jet can be absorbed 
by the accretion disk radiation field to create electron-positron 
pairs. As IC scattering, this process is also highly anisotropic. 
We estimate the opacity for photons emitted along the jet axis as 
a function of energy and z following Becker & Kafatos (1995). 
For this calculation, we also approximate the disk photon field 
at fixed radius as monoenergetic. 

5.2. Parameters of the fit 

In Table [3] we list the values of the model parameters that we 
keep fixed during the fit. Some of them, such as the black hole 
mass and the distance to Earth, are inferred from observations; 
we take their values from the literature about the source. For 
the remaining parameters, we have no information about XTE 
Jill 8+480 (such as the bulk Lorentz factor), or they are spe- 
cific to our model. For them, we adopt typical values for other 
microquasars, or simply estimates. 

The parameters in Table [4] were varied to obtain the best- fit 
models. The inner radius /?,„ and the temperature T max deter- 
mine the spectrum of the accretion disk. The value of g accr fol- 
lows from Eqs. ([TJ and (l44l . The rest of the parameters of the 
jet model (see Table [5) are calculated from the fixed and free 
parameters using the equations of Sect. 12.11 



7 Eqs. (15) and (16) of Bottcher et al. (1997) fail for values of their 
parameter a very close to unity. Since this interval does not add signifi- 
cantly to the value of the cooling rate, we disregard it in the calculation. 
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Table 3. Values of the fixed parameters of the model. 



Parameter 


Symbol 


Value 


Distance to Earth 


d 


1.72 kpc" 


Black hole mass 




8.5M S ° 


Disk outer radius 




7 X lOXrav* 


Jet viewing angle 


0jet 


30° 


Jet bulk Lorentz factor 


Tjet 


2 


Jet injection distance 


zo 


50R gmv 


Jet termination distance 


•£end 


10 12 cm 


Jet initial radius 


ro 


O.lzo 


Ratio L re l/Lj et 


q K \ 


0.1 


Magnetic field decay index 


m 


1.5 


Particle injection spectral index 


a 


1.5 


Acceleration efficiency 


1 


0.1 



Notes. (fl) Gelino et al. (2006). (i ° Chaty et al. (2003). 

Table 4. Values of the free parameters for the best-fit models. 



Parameter 


Symbol 


Value 






2000 


2005 


Disk inner radius 




164fl grav 


44tf grav 


Disk maximum tempera- 


T 

1 max 


22.4 eV 


46.5 eV 


ture 








Ratio Lj et / L accr 


<?jet 


0.16 


0.16 


Ratio t/s/f/k at Zacc 


P 


0.5 


0.85 


Ratio L p /L e 


a 


12.2 


25.5 


End of acceleration region 


Zend 


8.2 x 10 9 cm 


1.5 x 10 10 cm 


Minimum energy primary 


F 

^min 


86 ni( pe) c 2 


150 m (PiC )C 2 


protons and electrons 









Table 5. Best-fit model parameters calculated from the free pa- 
rameters. 



Parameter 


Symbol 


Value 






2000 


2005 


Accretion power 


^accr 


2.6 x 10 38 erg s 


-l 


9.9 x 10 37 erg s _I 


Jet power 




2.1 x 10 37 ergs 


-l 


8.1 x 10 36 erg s" 1 


Power relativistic 


Lp 


1.9 x 10 36 ergs 


-l 


7.8 x 10 35 erg s _I 


protons 

Power relativistic 


L e 


1.6 x 10 35 ergs 


-l 


3.1 x 10 34 erg s" 1 


electrons 










Magnetic field jet 


B 


1.3 x 10 7 G 




7.9 x 10 6 G 


base 










Base of the acce- 


Zacc 


2.8 x 10 s cm 




1.6 x 10 8 cm 


leration region 











5.3. Best-fit spectral energy distributions 

Figure[8]shows the best-fit SEDs for the 2000 and 2005 outbursts 
of XTE Jl 118+480. We obtain X l = 1-99 and X \ = 0.56 for 
the chi-squared per degree of freedom, respectively. The best-fit 
parameters are listed in Tables [3]|5j 

The value of the maximum temperature of the disk agrees 
with previous works (see for example McClintock et al. 2001, 
Markoff et al. 2001, Chaty et al. 2003, and Maitra et al. 2009). 
In the case of the 2000 outburst, T max is tightly constrained by 



the data, which clearly indicate the position of the peak of the 
multicolor black-body component. The value of the inner radius 
of the disk is not well-constrained by the observations, and de- 
pends on other details of the model as we discuss below. There 
is no agreement in the value of this parameter reported in the 
literature. 

We emphasize that our disk model is simple. We do not in- 
clude effects such as irradiation of the outer disk, or the transi- 
tion to an ADAF in the surroundings of the black hole. We do 
not attempt, therefore, to tightly constrain any characteristic pa- 
rameter of the disk, but just to account roughly for the thermal 
component observed in the SED. 

The typical accretion power of XRBs in low-hard state is 
L acc , ~ 0.01 - O.lLEdd- The value of g accr « 0.08 that we ob- 
tain from for the 2005 outburst is within this range. It is higher 
for the 2000 outburst, q accr ss 0.22, mainly because the disk in- 
ner radius is larger. Once r max and 9^ are fixed, R[ n determines 
the normalization of the spectrum. The maximum temperature 
of the disk for the 2000 outburst is well-constrained; the inclina- 
tion angle, however, is unknownQ The best-fit value of R m then 
depends on the value chosen for 9d, which we assumed to equal 
30°. Given the strong dependence of the accretion power on R m 
(see Eq. l44l i. small variations in this parameter yield significant 
changes in q. dccr . 

According to our modeling, there are no great differences in 
the physical conditions of the jets during the two outbursts. The 
radio and X-ray emission is fitted by the synchrotron spectrum 
of primary electrons, plus some contribution at low energies of 
secondary pairs created by yy annihilation. The IR-optical-UV 
range has significant contribution from the accretion disk. The 
synchrotron emission of secondary particles is negligible (except 
in the case mentioned above), as well as Bremsstrahlung radia- 
tion of primary electrons. The IC scattering off the jet photon 
field by primary electrons contributes in a narrow energy range 
at ~ 10 GeV in the case of the 2000 outburst. The SED above 
~ 1 GeV is completely dominated by gamma rays from the de- 
cay of neutral pions created in proton-proton collisions. 

The attenuation factor e~ T w as a function of energy and z is 
plotted in Fig. [9] The main source of absorbing photons is the 
accretion disk; the internal radiation field of the jet only adds 
a "bump" at high energies. The optical depth is large only near 
the base of the acceleration region. Gamma rays with energies 
10 GeV S £ y S 1 TeV are mostly absorbed in this zone. The de- 
velopment of electromagnetic cascades, however, is suppressed 
by the strong magnetic field (see Pellizza et al. 2010). 

The total luminosities are unmodified by absorption. As dis- 
cussed in Sect. 13.21 there are many high-energy protons that pro- 
duce gamma rays through proton-proton collisions outside the 
acceleration region. This radiation escapes unabsorbed because 
the density of disk photons is low at high z. 

The sensitivity curves of the Large Area Telescope (LAT) on 
board the Fermi satellite, the Major Atmospheric Gamma-ray 
Imaging Cherenkov Telescope II (MAGIC II), and the predicted 
for the future Cherenkov Telescope Array (CTA) are also plot- 
ted in Fig. [8] According to our results, a future outburst of the 
source with emission levels comparable to those of 2000 and 
2005, would be detectable in gamma rays by ground-based ob- 
servatories such as MAGIC II and CTA. 



8 There are estimates of the orbital inclination for XTE Jl 1 18+480, 
e.g. Gelino et al (2006). However, there is no compelling reason to as- 
sume that the disk lies in the orbital plane of the binary. See, for ex- 
ample, Maccarone (2002), Butt et al. (2003), and Romero & Orellana 
(2005) for a discussion of misaligned microquasars. 
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Fig. 8. Best-fit SEDs for the outbursts of 2000 (left) and 2005 (right) of XTE Jl 118+480. The various dotted curves labeled 
under the legend "synchr." correspond to the synchrotron radiation of primary electrons (the most luminous component), protons, 
and secondary particles (pions, muons, and electron-positron pairs). The dash-dotted curves labeled "IC" correspond to the SSC 
luminosity of primary electrons (red), and the external IC luminosities of primary electrons (violet) and muons (magenta). We also 
show the sensitivity curves of Fermi-LAT (1 yr exposure, 5<x), MAGIC II (50 h exposure), and CTA (50 h exposure). Figures are 
available in color in the electronic version of the manuscript. 
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Fig. 9. Total attenuation factor exp(-r rr ) caused by photon-photon annihilation for the best-fit SEDs of the 2000 (left) and 2005 
(right) outbursts of XTE Jl 1 18+480. Figure available in color in the electronic version of the manuscript. 



To date, no low-mass X-ray binaries have been detected with 
Cherenkov telescopes. Their observation is difficult because in 
general they are transient sources. Negative detections of four 
low-mass microquasars with HESS were reported in Chadwick 
et al. (2005). Three of them, however, were in the high-soft state; 
there are no available simultaneous X-ray data for the fourth 
source, GX 339-4, but it was apparently in a low-luminosity 
state. 

In the context of the model presented here, observations 
of low-mass microquasars at very-high energies would help 
us to constrain the hadronic content of the jets, since above 
~ 100 GeV the predicted emission is completely generated by 
proton-proton interactions. 

In addition, no low-mass X-ray binaries have been detected 
at high-energies. There is one Fermi gamma-ray source, 1FGL 
J1227. 9-4852, that might be the counterpart to the bright low- 



mass X-ray binary XSS J12270-4859, but the association re- 
mains unclear (Falanga et al. 2010, Hill et al. 201 1). 

The detectability in the Fermi-LAT band basically depends 
on the position of the synchrotron cutoff. We fix 77 = 0.1 for 
the acceleration efficiency, which yields a high maximum energy 
for the electrons, and the synchrotron emission extends into the 
MeV energy range. Observations with Fermi, then, could help 
us to investigate the efficiency of particle acceleration in jets of 
microquasars. 

6. Concluding remarks 

We have developed an inhomogeneous lepto-hadronic jet model, 
and presented some general results for the relativistic particle 
distributions, cooling rates, and the jet radiative spectrum. The 
model can satisfactorily reproduce the observed SED of the low- 
mass microquasar XTE Jl 1 18+480 in low-hard state. The best- 
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fit SEDs have revealed some interesting consequences of the in- 
jection of relativistic protons in an inhomogeneous extended re- 
gion. In particular, we have demonstrated that the gamma-ray 
emission can escape absorption, and that the source might be 
detectable when a new outburst of similar characteristics occurs. 

With some improvements, the same jet model can be applied 
to active galactic nuclei and gamma-ray bursts. The main modi- 
fications concern the calculation of the particle distributions. Eq. 
( TITT i would have to be modified to a covariant version, which 
could then be applied to outflows with very high bulk Lorentz 
factors. In addition, within environments where the radiation 
field is ultra-dense (as for the jets of long gamma-ray bursts), 
the coupling among the kinetic equations for particles and pho- 
tons could not be ignored. A further generalization of the trans- 
port equation would be to include time dependence, allowing 
the study of flares. Other interesting improvements to the model 
would be in the treatment of the jet kinematics and energetics, 
to develop a jet model that consistently couples the dynamical 
and radiative aspects. All these topics are part of our research 
project, and will be addressed in future works. 

Acknowledgements. G.S.V. thanks Iva Snidaric for the data on the sensitivity 
of MAGIC H. This work was supported by CONICET grant PIP 0078, by and 
the Spanish Ministerio de Ciencia e Innovacion (MINCINN) under grant AYA 
2010-21782-C03-01. 



References 

Abdo, A. A. et al. 2009, Science, 326, 840 
Albert, J., et al. 2007, ApJ, 665, L51 

Aharonian, F. A. 2004, Very High Energy Cosmic Gamma Radiation. World 

Scientific Publishing, Singapore 
Atoyan, A. M., & Dermer, C. D. 2003, ApJ, 586, 79 
Becker, P. A., & Kafatos, M. 1995, ApJ, 453, 83 
Begelman, M. C, Rudak, B., & Sikora, M. 1990, ApJ, 362, 38 
Blandford, R. D., & Konigl, A. 1979, ApJ, 232, 34 
Blandford, R. D., & Begelman, M. C. 1999, MNRAS, 303, L1-L5 
Blumenthal, G. R., & Gould, R. J. 1970, Rev. Mod. Phys., 42, 237 
Brocksopp et al. 2010, MNRAS, 404, 908 

Bosch-Ramon, V, Romero, G. E., & Paredes, J. M. 2006, A&A, 447, 263 

Bottcher, M., & Schlickeiser, R. 1997, A&A, 325, 866 

Bottcher, M., Mause, H., & Schlickeiser, R. 1997, A&A, 324, 395 

Butt, Y. M., Maccarone, T. J., & Prantzos, N. 2003, ApJ, 587, 748 

Chadwick, P. M., Dickinson, H. J., Latham et al. 2005, 29th ICRC, 4, 263-266 

Chaty, S. et al. 2003, MNRAS, 346, 689 

Chodorowski, M. J., Zdziarski, A. A., & Sikora, M. 1992, ApJ, 400, 181 

Dermer, C. D., & Schlickeiser, R. 1993, ApJ, 416, 458 

Dove, J. B., Wilms, J., & Begelman, M. C. 1997, ApJ, 487, 747 

Duro, R. et al. 2011, A&A, 533, L3 

Esin, A. A. et al. 2001, ApJ, 555, 483 

Falanga, M., de Martino, D., Bonnet-Bidaud, J. M. et al. 2010, Proceedings 
of the 8th INTEGRAL W orkshop, Dublin, Ireland. Published online at 
http://pos.sissa.it/cgi-bin/reader/conf.cgi?confid=l 15 id. 54 

Fender, R. P. et al. 2001, MNRAS, 322, L23-L27 

Frank, J., King, A. & Raine, D. 2002, Accretion Power in Astrophysics, 3rd 

Edition. Cambridge University Press, Cambridge 
Gaisser, T. K. 1990, Cosmic Rays and Particle Physics. Cambridge University 

Press, Cambridge 
Gallo, E., Fender, R., Kaiser, C. et al. 2005, Nature, 436, 819 
Gallo, E., Fender, R. P., & Pooley, G. G. 2003, MNRAS, 344, 60 
Gelino, D. M. et al. 2006, ApJ, 642, 438 

Ghisellini, G., Maraschi, L., & Treves, A. 1985, A&A, 146, 204 

Gould, R.J., & Schreder, G.P 1967, Phys. Rev., 155, 1404 

Heinz, S. & Grimm H. J. 2005, ApJ, 633, 384-391 

Hill, A. B., Szostek, A., Corbel, S. et al. 2011, MNRAS, 415, 235-243 

Hynes, R. I. et al. 2000, ApJ, 539, L37-L41 

Hynes, R. I. et al. 2006, ApJ, 651, 401 

Kaufman Bernado, M. M., Romero, G. E., & Mirabel, I.F. 2002, A&A, 385, 
L10-L13 

Kelner, S. R., Aharonian, F. A., & Bugayov, V. V. 2006, Phys. Rev. D, 74, 034018 
Kelner, S. R., & Aharonian, F. A. 2008, Phys. Rev. D, 78, 034013 
Khangulyan, D., Aharonian, F. A., & Bosch-Ramon, V. 2008, MNRAS, 383, 467 
Komissarov, S., Barkov, M., Vlahakis, N. & Konigl A. 2007, MNRAS, 380, 51 



Krolik, J. H. 1999, Active Galactic Nuclei: from the central black hole to the 

galactic environment, Princeton University Press, Princeton 
Laurent, P. et al. 201 1, Science, 332, 438 

Lipari, P., Lusignoli, M., & Meloni, D. 2007, Phys. Rev. D, 75, 123005 
Lyubarsky, Y. 2009, ApJ, 698, 2, 1570 
Maccarone, J. T. 2002, MNRAS, 336, 1371 
Maitra, D. et al. 2009, MNRAS, 398, 1638 

Malzac, J., Beloborodov, A. M., & Poutanen, J. 2001, MNRAS, 326, 417 
Markoff, S., Falcke, H., & Fender, R. 2001, A&A, 372, L25-L28 
McClintock, J. E. et al. 2001, ApJ, 555, 477 

McConnell, M. L., Zdziarski, A. A., Bennett, K. et al., 2002, ApJ, 572, 984 

Miller, J. M., et al. 2002, ApJ, 577, L15-L18 

Mirabel, I. F, & Rodriguez, L. F. 1999, ARA&A, 37, 409 

Miicke, A., Engel, R., Rachen, J.P., Protheroe, R.J. & Stanev, T. 2000, Comm. 

Phys. Comp., 124, 290 
Narayan, R., & Yi, I. 1994, ApJ, 428, L13 
Narayan, R., & Yi, I. 1995, ApJ, 444, 231 

Pellizza, L. J., Orellana, M., & Romero, G. E. 2010, IJMPD, 19, 671 
Pooley, G. G. 2005, Astron. Telegram, 385 

Poutanen, J., Krolik, J. H., & Ryde, F. 1997, MNRAS, 292, L21-L25 
Remillard, R., Morgan, E., Smith, D., & Smith E. 2000, IAU Circ, 7389, 2 
Remillard, R., Garcia, M., Torres, M. A. P., & Steeghs, D. 2005, Astron. 

Telegram, 384 
Reynoso, M. M., & Romero, G. E. 2009, A&A, 493, 1 
Reynoso, M. M., Medina, M. C, & Romero, G. E. 201 1, A&A, 531, A30 
Rieger, F. M., Bosch-Ramon, V., & Duffy, P. 2007, Ap&SS, 309, 119 
Romero, G. E. et al. 2003, A&A, 410, L1-L4 
Romero, G. E.,& Orellana, M. 2005, A&A, 439, 237 
Romero, G. E., Christiansen, H. R., & Orellana, M. 2005, ApJ, 632, 1093 
Romero, G. E., & Vila, G. S. 2008, A&A, 485, 623631 
Romero, G. E., Vieyro, F. L., & Vila, G. S. 2010, A&A, 519, A109 
Sabatini, S. et al. 2010 , ApJ Letters, 712, L10-L15 

Schlickeiser, R. 2002, Cosmic Ray Astrophysics. Astronomy and Astrophysics 

Library Series, Springer, Berlin 
Sell, P. H., Heinz, S., Calvelo, D. E., Tudose, V, Soleri, P. et al. 2010, ApJL, 

719, L194-L198 
Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337 

Spruit, H. C. 2010, The Jet Paradigm, Lecture Notes in Physics, Vol. 794, 

Springer- Verlag, Berlin Heidelberg 
Tavani, M. et al. 2009, Nature, 462, 620 

Tchekhovskoy, A., McKinney, J. C, & Narayan, R. 2008, MNRAS, 388, 551 

Tchekhovskoy, A., Narayan, R., & McKinney, J. C. 2010, New A, 15, 749 

Tomsick, J. A. et al. 2009, ApJ Letters, 707, L87-L91 

Torres, D. E, & Reimer, A. 201 1, A&A, 528, L2 

Vieyro, F. L., & Romero, G. E. 2012, in preparation 

Vila, G. S., & Romero, G. E. 2010, MNRAS, 403, 1457 

Yuan, E, Cui, W., & Narayan, R. 2005, ApJ, 620, 905 

Zurita, C, et al. 2005, Astron. Telegram, 383 

Zurita, C, et al. 2006, ApJ, 644, 432 



12 



